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Abstract 

New approach to the problem of multichannel continuum spectrum of three-cluster systems composed of an s- 
cluster and two neutrons is suggested based on the discrete representation of a complete basis of allowed states of 
the multiparticle harmonic oscillator. The structure of the eigenfunctions and behavior of the eigenvalues of the 
three-cluster norm kernel are analyzed. Classification of the eigenvalues of the three-cluster systems with the help 
of eigenvalues of the two-body subsystem is suggested. Asymptotic boundary conditions for a three-cluster wave 
function in the continuum consistent with the requirements of the Pauli principle are established. Such asymptotic 
behavior corresponds rather to subsequent decay of the three-cluster system than to the so-called "democratic decay" 
associated with the hyperspherical harmonics. The '^H-l-n -I- n configuration of the ®H nucleus is considered in detail. 
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1. Introduction 

During last years the superheavy hydrogen isotope has become an object of numerous experimental 
[112131415] and theoretical [6 7,8,9,10 researches. The overwhelming majority of these investigations aimed 
at finding the energy and width of the '^H resonance states. Unfortunately, even for the ground-state 
the experimental values for energy and width are significantly distinct and vary from Eif — 1.8 ± 0.1 MeV 
with r < 0.5 MeV ^ to Er ^ 5.5 ± 0.2 MeV with T = 5.4 ± 0.5 MeV [5;. Detailed discussion of the 
available experimental data is given in Ref. [llj with the conclusion that it is very difficult to find a non- 
contradictory scenario for all experimental data. Theoretical predictions are also different: in Refs. |8l9ll0j 
a resonance ,P — 1/2"'" is observed at ~ 1.5 MeV, while in Rcfs. [6|7I nearly twice that energy was 
obtained. As for the model, the most commonly used approaches treat the nucleus as a three-cluster 
system composed of a triton and two neutrons. Such assumption is justified by the fact that '^H is bound 
by 8.48 MeV [12], whereas ^n and ''H subsystems are known to be unbound. Mainly, different three-cluster 
models fall into two groups - macroscopic and microscopic. In macroscopic models clusters are considered to 
be structureless particles and cluster-cluster interactions are approximated by some local potentials which 
are fitted to reproduce relevant data on the cluster-cluster systems. As for the Pauli exclusion principle, it 
is usually simulated either with an additional repulsive potential between clusters or with orthogonalizing 
pseudopotential containing the operators of projection onto the forbidden states. The first approach was 
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used in Refs. [619] , where the blocking of the Pauh- forbidden states was reahzed by introducing a repulsive 
core in the s-wave ^H-neutron potentials. However, the choice of such a repulsive potential is quite ambiguous 
and within this approximation a complete and accurate exclusion of the forbidden states is not ensured. 
Furthermore, macroscopic three-cluster model rests on the assumption that cluster-cluster interactions are 
not affected by the presence of the third cluster, but this is by no means always the case. As regards the 
elimination of the Pauli- forbidden states with the help of orthogonalizing pseudopotential [13' , it could result 
in exclusion of some Pauli-allowed states. For example, in Ref. [T3] it was shown that some of the states of 
3a-system obtained by the above-mentioned method could be regarded as spurious, which they were not. 
Moreover, the latter states correspond to the most important shell-model configurations of the ^^C nucleus. 
Hence, the Pauli-allowed space should be carefully defined, in order not to exclude the dominant components 
for the realistic description of the three-cluster systems. For this purpose, one can define the three-cluster 
forbidden states referring to the harmonic oscillator wave functions of the microscopic two-cluster subsystem, 
as was suggested in Ref. [M]. This idea rested on the hypothesis that the antisymmetric three-cluster wave 
function should be orthogonal to the forbidden states of any two-cluster subsystem jl5] . This method seems to 
cope well with the task of accurate removal of the Pauli-forbidden states in a three-cluster system. However, 
it does not allow to distinguish the Pauli-allowed functions belonging to the same SU(3)-multiplet, i.e. the 
well-known problem of SU(3)-degeneracy of three-cluster states could not be resolved within such approach. 
Moreover, the elimination of the Pauli-forbidden states does not exhaust all exchange effects. Essential part of 
such effects is directly relevant to the eigenvalues of the antisymmetrization operator. The latter eigenvalues 
are not identical and determine the realization probability of the corresponding Pauli-allowed basis states in 
the wave function of the cluster system. Involvement of the eigenvalues of the antisymmetrization operator 
in the Schrodinger equation leads to changing the relative kinetic energy as clusters approach each other. 
Consequently, clusters are shown to experience an effective repulsion or attraction arising from the kinetic 
energy operator modified by the Pauli principle T6j. Such an effective interaction substantially affects the 
dynamics of the cluster-cluster interaction and can, on occasion, produce resonance behavior of the scattering 
phase or even a bound state in compound nuclear system [17 . As was shown in Refs. jl6|17j . the largest 
eigenvalues correspond to those basis states that are dominant in the discrete-spectrum states of the binary 
cluster system and in continuum states of small above-threshold energy. Also in Ref. [TS| it was observed 
that the probability of the presence of a cluster configuration in the lowest basis function for a binary cluster 
system is proportional to the eigenvalue of the isolated configuration. Hence, the eigenvalues of the Pauli- 
allowed states contain a wealth of information on a compound system composed of clusters. Meanwhile, 
the eigenvalues, along with corresponding Pauli-allowed states, depend only on the assumed internal cluster 
functions, not on the cluster-cluster potential, etc. Up to our knowledge, analysis of the structure of the 
eigenfunctions of the antisymmetrizer and behavior of its eigenvalues has never been performed for three- 
cluster systems, although we believe that it could help to establish some important laws that govern a 
three-cluster decay of nuclei composed of s-clusters. 

The main difficulty in studies of resonances in a three-cluster system consists in the formulation of the 
correct asymptotic boundary conditions for a wave function in the continuum. Such boundary conditions 
should ensure continuous transition from the region of small intercluster distances, where the Pauli exclusion 
principle is of first importance, to the asymptotic region where the scattering matrix elements are produced. 
Among currently available microscopic studies of ^H j7|8|10j , an explicit wave function representation of the 
scattering states have been employed only in Ref. [1^ . In Ref. [7] the method of analytic continuation in 
the coupling constant was used for study of resonance states in ^H, while in Ref. [8] the complex scaling 
method was applied for the same purpose. It should be noted that although authors of TU] performed 
calculations within the same model as we do, namely within an Algebraic Version of the resonating-group 
method (AVRGM) [19|, our approaches essentially differ in such aspects. To classify the three-cluster states, 
and enumerate the channels in the three-cluster continuum, hyperspherical harmonics were used in Ref. 
[lOj , while we define the Pauli-allowed harmonic oscillator basis functions in the Fock-Bargmann space and 
classify them with the use of the SU(3) symmetry indices. The latter seems more appropriate, because 
the second-order Casimir operator of the SU(3) group commutes with the operator of permutation of the 
nucleon position vectors. Hence, the SU(3) symmetry indices naturally appear as the quantum numbers of 
the eigenfunctions of the antisymmetrization operator. The hypermomentum, contrastingly, can not serve 
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as quantum number of Pauli-allowed basis functions, because the Pauli exclusion principle mixes the states 
with different values of hypermomentum. Hence, restriction for the maximum value of hypermomentum K 
included in calculation K < K^ax leads to spoiling the Pauli-allowed basis functions corresponding to the 
number of oscillator quanta v > K^ax, and this effect is enhanced with increasing i/. In support of this 
conjecture it was observed in Ref. [20j that an accurate description of the three-body asymptotics requires 
bases with large hypermomenta. 

In the present paper we have shown that correct asymptotic boundary conditions can be found employing a 
complete basis of the Pauli-allowed harmonic-oscillator states (classified with the use of the SU(3) symmetry 
indices and defined in the Fock-Bargmann space) along with their eigenvalues. Asymptotic behavior of basis 
functions consistent with the requirements of the Pauli principle gives an indication of possible decay channels 
of a three-cluster nucleus and allows us to specify the most important decay channels. Such asymptotic 
behavior correspond rather to subsequent decay of the three-cluster system than to the so-called " democratic 
decay" associated with the hyperspherical harmonics, which are widely used for the description of three- 
cluster systems (see Refs. [619110] ). 

In Section [2] the theoretical grounds of the AVRGM are outlined and the structure of the norm kernel 
for the number of three-cluster systems composed of s-cluster and two neutrons is discussed. Section [3] is 
devoted to the analysis of the eigenvalues and eigenfunctions of the norm kernel of the '^X = "^~^X + n + n 
systems {A < 6). The asymptotic behavior of the coefficients of the expansion of the three-cluster continuum 
wave function in the SU(3) basis is established in Section [4l In Section [5] the Pauli-allowed states of the 
^H^'^H+n + n system are considered in detail. The most important decay channels of the nucleus are 
specified and illustrated with the results of numerical calculations. Concluding remarks are given in Section 
|6l Definition and some properties of the Kravchuk polynomials are given in Appendix [A] 

2. Theoretical background and norm kernel 

Following the resonating-group method (RGM) [21], it will hereafter be supposed that the considered 
nuclear systems consist of three clusters. An RGM wave function is built in the form of an antisymmetrized 
product of cluster internal wave functions and a wave function of their relative motion. The internal wave 
functions of the clusters are and the wave function of relative motion of the clusters, which depends 

only on two Jacobi vectors of the considered three-cluster system, is found by solving an integro-differential 
equation. The latter is obtained by substitution of the RGM wave function into the Schrodinger equation 
followed by integration with respect to single-particle coordinates. The integro-differential equation can be 
transformed into the set of linear equations by expanding the wave function of the cluster relative motion 
into the compete basis of the Pauli-allowed harmonic oscillator states, as the AVRGM suggests. Another 
important simplification can be achieved by transformation from the coordinate space to the space of complex 
generator parameters (the Fock-Bargmann space [22i), in which basis functions arc of an especially simple 
form and are expressed via powers of complex vectors. Thus, the wave functions of the considered discrete 
representation are reduced to power series with an infinitely large convergence radius. The validity of this 
statement is indicated by the fact that all the wave functions in the Fock-Bargmann space are entire and 
analytic, and, therefore, series of these functions over powers of complex vectors converge in any finite region 
of the complex space. 

First and foremost, the AVRGM calls for the construction of the complete basis of the Pauli-allowed 
harmonic oscillator states and their classification. This is accomplished by solving the eigenvalue and eigen- 
function problem for the norm kernel, i.e., the overlap integral of the two Slater determinants composed of 
the single-particle orbitals: 



^ Here we shall assume the intrinsic cluster wave functions to be the simplest functions of a translation-invariant shell model 




(1) 
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Here integration is performed over all single-particles vectors, {Rj} (or {Sj}) identifies the set of three 
complex vectors determining position of the center-of-mass of clusters in the Fock-Bargmann space. For the 
spatial part of the single-particle wave functions we used the modified Bloch-Brink orbitals: 



<^(r») = ^ exp l^--rf + x/2(Rj • r,) - -Rjj , leA^, 

where Aj is the number of nucleons in the jth cluster. Each of these orbitals is an eigenfunction of the 
coordinate operator f^: 



which is defined in the Fock-Bargmann space and corresponds to the eigenvalue r^. and fjj are vectors 
of coordinate and momentum, respectively. At the same time, orbital 4>{vi) is the kernel of the integral 
transform from the coordinate representation to the Fock-Bargmann representation |i22j and the generating 
function for the harmonic-oscillator basis [55] . The center-of-mass motion are factored out (and dropped out 
from now on) by transition from generator parameters {Rj} to the Jacobi vectors: 

R-cm = —j= (^iRi + A2R.2 + ^sR-a) , 
V A 



^= MM^M)( Ml^^ + A^^A b=,/7HZ^(R.-R3). 

As the result, the overlap of two Slater determinants composed of the modified Bloch-Brink orbitals generates 
a complete basis of the Pauli-allowed harmonic-oscillator functions along with their eigenvalues: 



/({S,}, {R,}) ^ /(a, b; a, b) - ^ A„^'„(a, b)^'„(a, b). 



Functions are defined in the Fock-Bargmann space and orthonormalized with the Bargmann measure 

diiB = exp{-(aa)} exp{-(bb)} 

n stands for the set of quantum numbers of basis functions. Quantum numbers and the structure of functions 
^„ will be discussed in the next section. 

Examine next the explicit expressions for the norm kernels of those three-cluster systems which comprised 
of an s-cluster and two neutrons, such as ^Hc=*Hc+7i -I- 71, ^H='^H+n + n, '^n —"^ n + n + n, — n + n + n, 
subsequently referred to in this article as "^X = "^^^X + n + n, A < 6. The norm kernels for all these systems 
can be written in general form: 

^•S'(„„)=o = cxp(aa) cosh(bb) — cxp(a+a) cosh(b+b) — 

— exp(a_a) cosh(b_b) 4- exp{(aoa)}, (2) 

^S'(„„)=i = exp(aa) sinh(bb) — cxp(a+a) sinh(b+b) — 

— exp(a_a) sinh(b_b) + exp{(aoa)}. (3) 

Here S(^nn) designates the spin of the two-neutron subsystem. At that the so-called " T-tree" of the Jacobi 
coordinates is considered, i.e., vector b describes relative distance between two neutrons, while vector a 
indicates the distance from the remaining cluster to the center of mass of the two-neutron subsystem. 
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Note that /S(„„)=o is symmetric with respect to the permutation of two valence neutrons, while Is^„^y=l is 
antisymmetric as it must. 

Each term in Eqs. ^ and ([3]) corresponds to a certain permutation of identical nucleons. The first 
exponent is associated with the identity permutation, the second and the third exponents appeared as the 
result of permutation of a valence neutron and a neutron belonging to the s-cluster, while the last term 
relates to the simultaneous permutation of two pairs of identical neutrons: 

a± = (^l-0a±iVib, b± - ±i%/i a + ib, ao = (l-i)a. (4) 

Finally, the parameter t is positive and equal to the square of tangent of the angle of rotation from one 
Jacobi tree to another, t takes the value 3/2 for ^He nucleus, 5/3 for the ^H, t = 2 for tetraneutron and t = 3 
for three-neutron system. As can readily be observed, parameter t decreases with increasing the number of 
nucleons in s-cluster. 

The eigenvalues of the norm kernel do not depend on the choice of the Jacobi tree, whereas the structure 
of the eigenfunctions does. Since the Pauli-allowed basis states of the '^X — "^~^X + n + n systems take the 
simplest form in the " T-tree" , it precisely this tree is best suited to the construction and analysis of these 
states. 

3. Eigenvalues and eigenfunctions of the norm kernel of the — ^~^X + n + n systems 

Let us discuss now the set of quantum numbers of the Pauli-allowed states ^„(a, b). For the three-cluster 
systems considered here, n includes the number of oscillator quanta v, the indices (A, fi) of their SU(3) 
symmetry, the additional quantum number k if there are several differing (A, /i) multiplets, the orbital 
momentum L and its projection M and, if necessary, one more additional quantum number a^. The latter 
is needed to label the states with the same L in a given (A, fi) multiplet. The spin of the two-neutron 
subsystem 5'(„„) also is an integral of motion for as long as the spin-orbital interaction is switched off. It is 
well known that the diagonalization of the norm kernel requires the basis to be labeled with the quantum 
indices (A, /x) of irreducible representations of the SU(3) group [24] . The eigenvalues of the norm kernel 
depend on the total number of the oscillator quanta and {X,fi), and do not depend on the angular momenta 
of the basis states. We shall restrict our consideration to the states with zero orbital momentum L = and 
positive parity. Hence, the number of oscillator quanta should be even and equal 2;/, quantum numbers L 
and M will be dropped from now on. 

By definition, the states 5'„(a, b) are the eigenfunctions of the antisymmetrization operator A: 

In two-cluster systems the eigenvalues of antisymmetrizer limit to unity ai v —^ oo and the deviations from 
the unity are due to the Pauli exclusion principle. Contrastingly, eigenvalues of three-cluster systems tend 
to eigenvalues of a two-cluster subsystem with increasing the number of oscillator quanta v. 

.-2™ooV-4m,2m). 

Coincidence of the eigenvalues of the three-cluster norm kernels with those of the two-cluster subsystem 
norm kernels in the limit of the large number of quanta has been already pointed out in 26], where 3a 
and ^^0+2a systems have been considered. However, no relation with the additional quantum number k 
labeling the states with the same SU(3)-indices has been established. This subject will be pursued further 
in the next section. Now let us proceed with the Pauli-allowed basis functions ^(2^"4yj 2/^)^ 

In the Fock-Bargmann space the latter functions are the superpositions of the eigenfunctions of the 
second-order Casimir operator: 

*(2^"4m,2m) J^' ^ X! -^^2™-4M,2M)fc'^^2™-4l^,2M)(^''^)' 
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1/— /J.— 1 

^(2i/-V,2^)fcl^' "-I - ^(2zy-4Ai,2p)fc ^(21/- V,2ai)'^^' 

m— /i 

Functions ■'/'^2i/-4/j 2/^) ^-"^^ ^^2"-4^^^) correspond to the same total number of oscillator quanta 

2v, but differ in the number of quanta along vectors a and b. That is the reason why the three-cluster 
states are always degenerate, even though only s-clusters involved, and the degree of the SU(3) degeneracy 
increases drastically with increasing the number of quanta. It was found in Ref. |25| that the eigenfunctions 
of the Casimir operator with powers of vectors a and b being fixed are expressible in terms of hypergeometric 
functions 2-^1(0!, /?; 7; Z), with the variable 

fabl2 



Z = 



a2b2' 

Namely, 

^?2::f,,2.)(a,b) ^ <:,t2,)[ab]^-a-^---^b--- X 

X2i^i (^v + ni + fi, — m + + 2fi+ -■,Z ] , 

4:=4r2!.)(-'b) =<:4tt;,)[ab]^-(ab)a-----^b--- x 

X2-F'i ^— + m + + 1, -rn + /i; -i^ + 2fi + -; Z 

Here -^(^2™-4m 2/i) ^^"^ ^{■iv-i'/ttli) normalization coefficients of the corresponding basis functions. 

The norm kernel Is(„„-, can be written down as a sum of partial norm kernels with definite values of the 
number of quanta 2i': 



(1 



In the Fock-Bargmann representation, at a given v the norm kernel ^ is always representable in the 



form of a sum of S'C/(3)-projected norm kernels l's'^^\- We shall deal with the relevant part of the norm 

(21^-4/^,2/^) 



kernel, />, ' and write it as 



4(1^=0^^'' = XI XI ^2^' 2™ ~ ^^l^'^' 2m - 2^ + 2m) x 

m—fi rh—Q 

I I 2m — 2/i 72m — 2fi-\-2rh , ( 2m — 2/^+2m 72m — 2fi \ 
^ \^(2i^-4:fi,2fi)'^(2iy-4fi.2fi) ^ ^(2i/-4m,2m) ^(21.-4^:2^)^ ' 

i^ — fi—l v — m — fi—l 

Is^Z^=l^^ = X X 2m + 1 - 2^|2;/, 2m + 1 - 2/i + 2m) x 

f j2m+l-2fj. T2m+l-2/j+2m . , 2m+l-2/j+2m T2m+l-2/j \ 
^ \^^(2iy-4^,2p)'V^(2iy-4^,2^) +''^(21^-4^,2/*) ''^(2iy-4^,2^)^ ' 

where 

(2i/, 2m - 2^j.\2iy, 2m - 2^ + 2m) = ^^.o (l + (5m,o'5/.,o(l - tf") - 



(t_l\ "'^ _ / (2i^-2m-2At)!(2m + 2m- 

\t-2j ^ ' y (2j/- 2m- 2m- 2^)!(2m- 

X2F1 ( -2i^ + 2m + 2/i + 2m, -2m + 2/i; 2m + 1; 



2~t 



{2v, 2m + 1 - 2^\2v, 2m + 1 - 2^ + 2m) = ^.^,0 + 



2t^ ^t-lV''.,_2)2(.-™-A)-i / (2^-2m-2/i-l)!(2m + 2m+l-2A.)! 

^ ' Y (2i/- 2?7i- 2m- 1 - 2^)!(2m+ 1 - 

X2F1 ( -~2v + 2m + 2/i + 1 + 2771, -2m -1 + 2/^; 2?7i + 1; 



Obviously, the Pauli- allowed basis functions ^'(2^"4^ 2^l)^, '"^"^ found by the diagonalization of the norm 
kernels /(.■|''~^'^'^^) at a given u. There are v ~ 2^ + 1 different states possessing SU(3) symmetry {2v — 
4/i, 2/x), but some of them correspond to zero eigenvalues and thus they are forbidden by the Pauli principle. 
The number of the latter functions depends on the system considered. Hence, solving the eigenvalue and 
eigenfunction problem for the norm kernel /(a, b; a, b) in the Fock-Bargmann space, we arrive at all the 
Pauli-allowed and the Pauli-forbidden functions of the "^X = "^~^X + n + n systems given in terms of 
orthogonal polynomials of a discrete variable. 

^A-2jr_|_„ , „ _ , , . A-1 ■ - 

'-(2iy-4M,2/i)fc 

eigenvectors ^(21^-4^1, 2/^)^ take simple analytical form 



As the eigenvalues A/2^_. „''^^"^" of the norm kernel approach limit values Aj, ^ corresponding 



*(2i/-4M,2/i);.^ 



(»' ^) ^ 71 (^'^2.-4^,2^) (ai, bl) + ^f2,/-4p.2M) (^2, bj)) , (5) 



^S^-M). ^ 7f ('^('2.-V,2m) (^1' '^l) - ^f2.-V,2M) (^2, b2)) . (6) 

Such asymptotical behavior takes place with the proviso that v ^ k, but as will be exemplified with 
system, the fulfilment of condition v > k + 5 would be ample. Note that k is the natural integer that falls 
in the range 2/1 < k < v. 

Here ai,bi and a2,b2 are the Jacobi vectors of the so-called "Y-tree", with vector bi (b2) describing 
relative distance between one of the valence neutron and "^~2x-cluster. Correspondingly, vector ai (a2) 
determines position of the remaining neutron relative to the center-of-mass of the '^"^X subsystem. Jacobi 
vectors of the "T-tree" are related to those of the "Y-tree" via unitary transformation: 

a = cosQ;ai_2 4- sinabi_2;b — =psinaai^2 i cosabi^2- 

Naturally, function ^^2i/-4Ai'^2/i)fc symmetric with respect to the interchange of the vectors (ai,bi) and 

S —1 

(a2, b2), while ^(-2^'"4p 2/^)4. antisymmetric, in so far as such operation corresponds to the permutation of 
the valence neutrons. 

The remarkable feature of asymptotical relations ([5]) and ^ lies in the fact that in the limit v ^ k 
the expansion coefficients ^^^Z4^ 2/^)^ ' ^^2i!'-4;i^Ji)fc '^^'^ identified with the Kravchuk polynomials of a 
discrete variable \T7\: 



DlZ-X%. ^^i''-^2,(2m -2^ + 1) 



dk-2ii 

2m-2M+l ^,/^ir(P) 10^ O., I i^\/P2m-2A.+ l ^ _ .^2 



p — sm a. 



The Kravchuk polynomials ^'^-2^^''^ ^ -^A*) ^'^^ specified on the interval 2^ < m < 2v — 2/i and orthogonal 
with weighting function Pm-2p, and norm dk-2fj.- An explicit expression for the Kravchuk polynomials, as 
well as some of their basic properties, are given in Appendix [X] Here it is worth noting that the Kravchuk 
polynomials are the discrete analog of the Hermitian polynomials. 
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So, the degree of the Kravchuk polynomials serves as the additional quantum number of the SU(3) 
degenerate three-cluster states. The Pauli-allowed basis states ^'fj""', ^ % can be arranged into branches 
and families, with all the states of a particular branch having common symmetry index /i and overlapping 
generously with corresponding asymptotic function ([5]) or ([6]), but differing in value of the first index A. 
The eigenvalues belonging to a given branch tend to the same limit value A^, with the number of quanta 
increasing. The branches which share limit eigenvalues are combined in the family of the eigenstates, which 
thus is completely determined by the degree k of the corresponding Kravchuk polynomial. 

Evidence for the importance of a particular family of the Pauli-allowed states can be found by analyzing 
the behavior of their eigenvalues with the increasing the number of oscillator quanta v. In the Introduction 
mention was made of the fact that the change in the kinetic energy of relative motion under the effect of 
the Pauli exclusion principle leads to an effective interaction of the nuclei as they approach each other. As is 
shown in Refs. |16I17| . the character of such interaction in two-cluster systems is determined by the behavior 
of the eigenvalues of the antisymmetrization operator. They are nonnegative, since they are proportional to 
the probability of the realization of the corresponding allowed basis state, and tend to unity with increasing 
distance between the nuclei involved. In the case where the eigenvalues of the antisymmetrization operator 
tend to unity from below, the effective interaction proves to be repulsive, but, if these eigenvalues tend 
to unity from above, one can say that there arises attraction induced by the Pauli exclusion principle. As 
indicated earlier, the eigenvalues A^2^~^^ 2^)t"^" °^ three-cluster norm kernel approach the eigenvalues 
of the two-cluster subsystem Aj, ^ with increasing the number of quanta ly. It would appear reasonable 
that families of the states with A(2j/-4p,2;^)fc < Afc would be suppressed at small which can naturally be 
considered as a demonstration of the action of effective forces of repulsion for small intercluster distances. 
Contrastingly, those families which characterized by A(2i^-4/j,2/x)fc > should become more favorable at 
small i^, which can be considered as effective attraction. The validity of these assumptions for the 
nucleus will be demonstrated in Section [5l 

It is worth noting that the effective interaction caused by the change in the kinetic energy of the relative 
motion of clusters under the effect of the antisymmetrization operator arises only between objects whose 
internal energy may change when they approach each otheiT^. Therefore, no such interaction would appear 
between two point-like nucleons. By this reason, the Pauli-allowed states in ■^n=n+n+n system will most 
likely have the asymptotic behavior different from 

4. Asymptotic equations for the expansion coefficients 

We seek the wave function of the considered three-cluster system in the form of expansion over the SU(3) 
basis of the Pauli-allowed states 

T,(B)(a,b)=^V^C:(^)vl/„(a,b). (7) 

n 

Expansion coefficients both of the discrete eigenstates with energy E,^ = — < and of the continuum 
eigenstates {C„(i?)} with energy E > are found by solving a set of linear equations 

Y,{n\H\n')Cn' - EAnCn = 0. (8) 

n' 

In Ref. [28| we have shown that in order to set the asymptotic boundary conditions for the expansion 
coefficients of a two-cluster wave function in the SU(3) basis, a basis with a different set of quantum numbers 
(the angular-momentum coupled basis) is required. The states of the latter basis are labeled with angular 
momenta of the clusters and of their relative motion. The equations ([8]) for the expansion coefhcients in the 
SU(3) basis remain coupled even in the asymptotic region, whereas the set of corresponding equations in the 
angular- momentum coupled basis is uncoupled at a large number of excitation quanta. The transformation 

^ Here, we imply the interaction of nuclei consisting of nucleons. 



8 



between the two is defined through a unitary matrix and can be found with the use of the integration 
technique in the Fock-Bargmann space developed by the present authors. However, the norm kernel of 
binary cluster system has a diagonal form only in the representation of the basis of its eigenfunctions, i.e., 
in the representation of the SU(3)-basis. A unitary transform of the basis disrupts the diagonal form of the 
norm kernel due to the difference between the eigenvalues for different SU(3) representations. This would not 
be the case if all the eigenvalues were equal. The nature of this breaking is that, unlike the functions of the 
SU(3) basis, those of the angular-momentum coupled basis are not eigenfunctions of the antisymmetrization 
operator and, therefore, are not invariant with respect to a permutation of the nucleons. The permutation 
mixes the angular-momentum coupled basis functions with the same number of quanta. However, as the 
latter increases, the degree of mixing decreases exponentially, and at large values of the number of oscillator 
quanta the norm kernel becomes practically diagonal in the angular-momentum coupled basis as well. Hence, 
in this region the asymptotic behavior of the expansion coefficients in the angular-momentum coupled basis 
can be defined (they are expressed in terms of Hankel functions of the first and second kind, and the 
scattering S'-matrix elements) and related to that of the expansion coefficients in the SU(3) basis. 

Contrastingly to the eigenvalues of the two-cluster norm kernels, those of the three-cluster norm kernel 
remain distinct even in the limit ly ^ oo. Existence of different limit eigenvalues Afc reflects the possibility for 
two of three clusters to be close to each other and far apart from the third cluster. Owing to this the unitary 
transformation from the SU(3)-basis to the other one is rendered possible only within a particular family of 
the Pauli-allowed states. Any transformation, which involves the states belonging to different families (for 
example, transformation to the three-cluster hypcrspherical harmonics), would disrupt the diagonal form of 
the norm kernel and, hence, is inappropriate to the occasion. 

Turning back to the asymptotic SU(3) basis functions jS]),!!]), it might be well to point out that the latter 
functions have a simple physical meaning. They reproduce relative motion of a two-cluster subsystem "^~^X 
occurring in a ground or an excited harmonic-oscillator state and a remaining neutron. Such asymptotic 
behavior gives an indication of possible decay channels of a three-cluster nucleus and allows us to specify 
the most important decay channels of the nucleus under consideration. 

In the limit v ^ k instead of asymptotical SU(3) basis functions 

^ (v^(V.'-4m,2m) (ai ' bi ) ± ^\2^!^i^,.2^,) (a2 , b2 ; 
it is advantageous to use the angular momentum basis functions 
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(ai,bi) ± (t>2u,ki.^2M)) 



defined as follows: 

[1/2] 

.2.,,(a,b)=7V^,,a^-^-'b^-'(ab)'.^ 



[=0 



2' ^n(;-0!(/-20!(ab)2' 



The states of this basis (referred to as "/-basis" in what follows) are labeled by the number of quanta 2i^, 
and the angular momenta / of the "^~^X subsystem coinciding with the angular momentum of the relative 
motion of this subsystem and a remaining neutrorf^. ^ is the normalization coefficient. 

For ly ^ k variables in ([8]) are separated and the latter set of equations is representable in the form of a 
sum of two sets of equations. One of them describes the relative motion of '^^^X subsystem and a neutron: 

-\V{^'^~k-l){2:.-k + l + l)Cl_,^, + U2iy~k+^^2E] Cl^,- 



- -^{2iy -k-l + 2){2v -k + l + 3)C^+i^fe ^ 0, (9) 
while the other one characterizes the "^^^X =^^2 x + n subsystem itself: 



^ Recall that only the states with total angular momentum L = are considered in this paper. 
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-J^-y{k-l + 2)ik + l + 3)C:+i,,+2 = 0. (10) 

We here put the nucleon mass, the Planck's constant and oscillator length equal to 1 for the sake of brevity. 

Thus asymptotically the three-cluster Schrodinger equation ([5]) can be reduced to a two-body-like mul- 
tichannel problem. Total energy E is distributed between energy e of the two-cluster subsystem "*~^X and 
energy E of the relative motion of the above-mentioned subsystem and a remaining neutron so that E+e = E. 
Note that Eq. ^ corresponds to the free motion of a neutron and "^"^X subsystem, while the matrix of Eq. 
^ contains the limit eigenvalues Xk of the allowed states and, for this reason, is not identical to the matrix 
corresponding to the kinetic energy operator of free motion of the '^~^X cluster and a neutron. In other 
words, it contains some effective cluster-cluster interaction derived from the kinetic-energy operator modified 
by the Pauli principle. It is precisely this interaction that causes the ^X system to decay via an intermediate 
stage, i.e.,'^X X + n -^^"^ x + n + n. Alternative decay channel "^X X -f ^ n ^^^^ X + n + n 

is much less favorable, because in this case Eq. (fTO)) would describe just a free motion of two neutrons and, 
thus, a sufficiently high energy is needed to make this decay channel open. The possibility for such decay 
channel to exist in is discussed in Section [5l 

From the above discussion it appears that the greater is the number of the families of the Pauli-allowed 
states invoked, the better is the description of the '^~^X subsystem and the less is its energy e. Failure to 
take into account a sufficient number of families results in too localized "^~^X subsystem and, hence, in too 
high energy s. Noteworthy also is the fact that in the limit of the large number of oscillator quanta v the 
kinetic energy operator does not couple the states belonging to the families corresponding to even values of 
quantum number k with those characterized by odd-numbered k. Such phenomenon takes place because even 
and odd values of k correspond to the states of the "^~^X subsystem with different parity, which becomes 
an integral of motion at large distance between the '^~^X subsystem and a neutron. 

Of course, the cluster-cluster interaction generated by the nucleon-nucleon potential should also participate 
in Eq. ([9|) providing for realistic description of the "^~^X subsystem. However, inclusion of the nucleon- 
nucleon interaction between nucleons belonging to different clusters could not change the structure of the 
Pauli-allowed states and in a qualitative sense affect asymptotic behavior of the expansion coefficients of the 
wave function over these states. Study of the effects related to the potential energy (with its exchange part) 
in the microscopic analytical approach described here is under its way, and the results will be published in 
a separate paper. 

Difference equations ^ and pO|) become Bessel equations in the limit of large v and k: 



(11) 



(^ + -:^ + ^^-^^^^)cUyk) = 0; yk^V2kTi. (12) 
\dyi ykdyk ^yi J 



Therefore, the asymptotic form of the expansion coefficients in the Z-basis can be conveniently written in 



terms of the Hankel functions Hj^-^,^ and the scattering S'-matrix elements. If the incoming wave is in the 



channel characterized by angular momentum I and energy E — Si oi the relative motion of ^X subsystem 
and a neutron, the expansion coefficients in this channel satisfy the asymptotic relation 



£ S,,H^+l, ( ^2{E e,) V4Z. - 2fc + 3 j iJ^/^ (v^V2fcT3) 
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Here index i enumerates different decay channels and its maximum possible value is equal to the number of 
different families considered. We emphasize that energy levels Si of the '^"^X subsystem do not necessarily 
belong to the discrete states, but can be continuum states as well. Even in the latter case it is possible to 
indicate the most favorable decay channels. For example, the "*~^X subsystem may not have any bound 
states, but could have a low-energy resonance. In such a case the "^^^X + n + n three-cluster system would 
most likely decay via resonance state of the '^~^X subsystem. 



5. The Pauli-allowed states of the ^H+n + n system 



Let us use a three-cluster configuration ^H-|-n -I- n of the nucleus with — 1/2+ to illustrate the 
validity of our conclusions. We shall focuss our attention on the states with positive parity, zero total 
orbital momentum and zero spin of the two-neutron subsystem. For the sake of brevity the above-mentioned 
quantum numbers will be dropped in what follows. 

In this case the norm kernel ([2]) takes the form: 



I = cosh(aa) cosh(bb) — cosh(a+a) cosh(b+b) 
— cosh(a„a) cosh(b_b) + cosh{(aoa)}. 



(13) 



where a±,b± and ao are defined by Eq. (jH) with t — 5/3. 

The norm kernel (fT3)) have to be projected to the states with definite indices of the SU(3) symmetry 
(2i/— 4/i, 2/i) and expanded over a complete set of the eigenstates ^'(2i/-4p,2/^)fc (^-i b) of the antisymmetrization 
operator: 



^ X! ^(2i^-V,2,j)fc*(2i/-4M,2Ai)fc- 
fe=2/i 



Eigenvalues of the norm kernel can be arranged into the branches and families quite as discussed in Section 
|3l Eigenvalues belonging to the first five families are given in Table [TJ 



Table 1 



Eigenvalues A(2i^— 4^,2^);. of the norm kernel for the system '^H-l-n + n 



k -- 
V fi - 



1 k = 2 

: At = ^i = 



1 



k = 3 
= fi = 



1 Pi- 



k -- 
fi - 



: 4 

:1 fi-. 



■2 pL- 



k ■- 
fi 



■■ 5 

:1 



1 1.8889 

2 1.4780 0.9541 0.7778 

3 1.3487 0.8814 0.8733 1.0895 1.0649 

4 1.3344 0.8862 0.8879 1.0636 1.0434 1.0161 0.9837 0.9753 

5 1.3334 0.8885 0.8888 1.0431 1.0377 0.9841 0.9847 0.9859 1.0139 1.0034 1.0072 

6 1.3333 0.8888 0.8888 1.0378 1.0371 0.9861 0.9870 0.9875 1.0088 1.0055 1.0048 

7 1.3333 0.8889 0.8889 1.0371 1.0370 0.9873 0.9876 0.9876 1.0056 1.0045 1.0042 

8 1.3333 0.8889 0.8889 1.0370 1.0370 0.9876 0.9876 0.9876 1.0044 1.0042 1.0041 

9 1.3333 0.8889 0.8889 1.0370 1.0370 0.9876 0.9876 0.9876 1.0042 1.0041 1.0041 

10 1.3333 0.8889 0.8889 1.0370 1.0370 0.9876 0.9876 0.9876 1.0041 1.0041 1.0041 



One can readily see from Table [T] that the norm kernel eigenvalues belonging to the kth family tend 
to the eigenvalues of a two-cluster subsystem with increasing the number of oscillator quanta ly; 



A 



'H=^H-(-n 
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Obviously, index k makes sense of the number of oscillator quanta accounted for by the subsystem. Of 
special note is the fact that Aj.^2fc+i > ^^ while Aj,^2fc < 1- This in itself is indicative of attraction between 
•^H and a neutron in the states with odd number of quanta k — 2k + 1 and repulsion in the states with 
even number of quanta k = 2k. Such a conclusion can be easily verified by analyzing the asymptotic Eq. 
(fTO)) . Moreover, the eigenvalues A(2y_4^_2p),j_|_j approach limit eigenvalues X2k+i fro™ above, contrastingly 
to A(2!y-4;2.2^),^ approaching X^j, from below. As was concluded in Section [Sj there are strong grounds for 
believing that families of the states characterized by odd values of quantum number k dominates in the 
wave function of the '"^H system. Table [2] illustrates the fact that the Pauli-allowed states 5'(2iy-4;^.2^t)fc ^^^^ 

Table 2 

Overlap integrals J 'I'(2i'— 4/j,2/j)fc ^('2y_4p 2fj) '^f^B versus the number of quanta i/ for the first five families of the Pauli-allowed 

states of the ^H. 

fc = l fc = 2 A: = 3 fc = 4 k = 5 

V fi = fj, = fJ. = 1 = /t = l ^ = /i = l fi = 2 fj, = At = l fi = 2 

1 0.9701 

2 0.9460 0.9924 1 

3 0.9893 0.9979 0.9934 0.8076 0.9915 

4 0.9996 0.9998 0.9992 0.9007 0.9886 0.9214 0.9629 1 

5 1.0000 1.0000 1.0000 0.9736 0.9983 0.9706 0.9900 0.9934 0.5887 0.9052 0.9915 

6 1.0000 1.0000 1.0000 0.9981 0.9999 0.9962 0.9979 0.9992 0.7299 0.9907 0.9886 

7 1.0000 1.0000 1.0000 0.9999 1.0000 0.9995 0.9998 1.0000 0.9294 0.9974 0.9983 

8 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9932 0.9996 0.9999 

9 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.9996 1.0000 1.0000 

10 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 

asymptotic form 

as A(2iy_4^_2/i)fc approach their limit eigenvalues. In Table [2] overlap integrals of the exact Pauli-allowed states 
with their asymptotic expressions ^^2i/-4^ 2/j)fc depending on the number of quanta 21/ are given for the first 
five families of the norm kernel eigenfunctions. As long as < 5, the states listed in Table [5] exhaust all 
possible basis functions allowed by the Pauli principle. As ly increases, new families of states characterized 
by fc > 5 also appear. However, eigenstates and eigenvalues belonging to such families are governed by the 
same laws as those tabulated in Tabled Furthermore, states of "odd" families, i.e., ^'(21^-4^,2/^)2^^^ appear 
to dominate over the states 4'(2i/-4/i.2/i)2^ 7 with small k prevailing. The first family 5'(2i/,o)fc=i pl^-ys a leading 
part in the wave function of the system, followed by the famihes ^'(2iy-4/i, 2^1)^=3 and 'i'(2u-4fi,2fj.)k=5- Such 
hierarchy among different families of the Pauli-allowed states of the is directly related to the behavior of 
the corresponding eigenvalues ^(21^-4^^2(1)2^^^ discussed above. An exact treatment of the antisymmetrization 
effects related to the kinetic energy exclusively was shown to result in an effective attraction of the clusters 
in those branches whose eigenvalues exceed their limiting values. 

From Table [1] and Table [5] we can observe that a complete space of the Pauli-allowed basis functions 
can be divided into "internal" region, where — k < 5, and asymptotic region {v — k > 5). In the latter 
region functions ^(2iy-4/i,2/j)fc take simple analytical form '^^2i'-4,(i 2^1)^ ^^'^ correspond to the decay of the 
nucleus into the subsystem and a neutron, with the cluster occurring in the state with k oscillator 
quanta. In the internal region the Pauli-allowed states have more complicated form, which is indicative of 
exchange effects involving all three clusters. 

The states of the fcth family have one more important feature, which makes small k predominating. 
The less is fc, the more smooth basis functions ^'(21^-4/^,2/1)^ are. Let us demonstrate this with the functions 
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belonging to the first family consisting of a single branch 5'(2i/-4/^,2/j)fc=i ■ The latter branch has a remarkably 
simple asymptotic form: 

"^tko). - ^ ■ ^(^^ {af~^(aibi) + ar-^(a2b2)} , 

and the exact function ^'(2;y-4/i, 2/^)^=1 becomes indistinguishable from the asymptotic one at = 5. It is 
easily comprehended that function '^^2i> 0)1 corresponds to the decay of the nucleus into the "^H subsystem 
being in the lowest oscillator shell model state and a neutron. Indeed, ^'(Ij^o)^ contains only the first power 
of vector bi or b2 evidencing that cluster and a neutron remain at the minimal distance from each other, 
which is consistent with the requirements of the Pauli exclusion principle. In turn, eigenvalues of the first 
family A(2,y,o)i tend to the lowest eigenvalue of the ''H norm kernel \k=i — 4/3. It is noteworthy also that 
the orbital angular momentum I of the relative motion of the neutron cluster and subsystem, as well as 
the orbital momentum of the itself is Z = 1 for the states of the first branch. 

As noted in Section [31 at > fc + 5 the expansion coefficients D'^^Z.^'^ 2^)1- ^^"^ Pauli-allowed basis 
functions ^'(21^-4/^,2/^)1. ^an be identified with the Kravchuk polynomials of a discrete variable 2m — 2/i: 



Dependence of the expansion coefficients of the first branch ^'(2iy~4/j, 2/^)^=1 on the number of quanta v is 
shown in Fig. [1] At a given v there are exists I' + l coefficients -D^^ 0)^-1 ■ expected, the latter coefficients 




Fig. 1. Expansion coefficients of the first branch versus m at different values of the number of quanta u. Curves: 1 

have one node at all v, because their asymptotic behavior is determined by the Kravchuk polynomials of 
the first order. Note that with increase in f coefficients D^^Z o)k-i become Hermitian polynomials of the first 
order as follows from (jA.ip . Evidently, the less is the number of nodes, the lower is the kinetic energy and 
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potential energy corresponding to a given branch. Therefore, famihes of the aUowed states associated with 
the low-order Kravchuk polynomials are more favorable. 

It should be pointed out that amongst eigenfunctions of the antisymmetrization operator of the system 
belonging to SU(3) representation (2i/, 0) there is a nodeless branch of eigenfunctions ^'(2;y,o)fe=o- However, 
it falls into the category of the Pauli-forbidden states and corresponds to zero eigenvalue. Functions of this 
branch can be written in analytical form for a given v: 

where N(^2v,o)^_ is the normalization coefficient, which tends to l/-\/2 with ly increasing. The expansion 
coefficients -D^;^o)fc-o given by the following expression: 



D 



2m 




{2iyy. /5 



The structure of '^{2u.o)k=o indicates that one of the valence neutrons occupies an s-state in the ''H subsystem, 
which is forbidden by the Pauli principle. 

Let us discuss now possibility for another decay channel of the to exist. At first glance it would 
seem that ^H^^ n+^H^ n + n+^ii decay should occur along with the ^H— > n+'^H^ n + n+^H decay. 
However, the structure of the Pauli-allowed functions is not compatible with such assumption. Basic functions 
corresponding to the "dineutron" decay of the ^H, if any, should look like 

^'(2i.-v,2p)2™(a,b) -> V^2T4il,2M)(^''=*)' M<w<i^-/i. (14) 
with eigenvalues 

^2iy-4A«, 2^1)2™ ~^ ^2m 



lim A(r>,,-A,, 9„~io™ Ao" — 1. 



1^ — 2/^ — >oc 

Asymptotic behavior of the simplest function relevant to the decay of the into the dineutron being in 
the lowest shell-model state and triton is of the form 

From this it readily follows that 
lim DI^^q) Synfl- 



In Table [3] eigenvalues of those basis states which most closely resemble functions are listed along 
with the weight of the state "0(^170) (a, b). 

Table 3 

Eigenvalues of the candidate for the "dineutron decay" function ^{2iy,o)„-o '-''^ ^'^'^ overlap 

^T^%) = / "^fL,0)„=o*(2i',0)„=o'^/^s versus number of quanta v. 



V 123456789 10 

A{2^,0)„=o 1-8889 1.4781 1.0895 1.0636 1.0139 1.0087 1.0027 1.0015 1.0008 1.0003 
■^(2^0) 0-8575 0.7304 0.8407 0.7327 0.8564 0.8510 0.8141 0.9518 0.8237 0.9406 



As may be inferred from this table, the dependence of the coefficient -C'j^~g-| from the number of quanta v 
is not monotone and even at = 10 the weight of the "dineutron decay" function is only 88%. Note that the 
weight of the simplest ^H+n function in the states of the first family runs up to 100% '&i v — h. As regards 
the branches 5*^21/ 0)2 >o ' ^'^i*^'^ describe the decay of the into the dineutron being in the state with 
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2m quanta and a triton, none of the Pauli-allowed states of the has such asymptotic behavior. Figure [5] 
depicts the expansion coefHcients D'^2i'0) -o branch, which is a candidate for the "dineutron decay" 

mode, at different values of the number of quanta ly. As evident from Fig. [21 coefficients -D?^ _^ have 




Fig. 2. Expansion coefficients ^ of the "dineutron decay" branch versus m at different values of the number of quanta 

u. Curves: 1 - ly = 5; 2 - u = 10; 3 -'u = 15. 

noticeable tail even at v = 15, indicating that function ^'(2iy,o)m=o ^^^^^ does not coincide with its asymptotic 
expression 5*^21/ o) -o ■ Furthermore, it has too many nodes and thus quite high energy is needed to excite 
such a mode. All the foregoing counts in favour of the conclusion that ^H— ^■^ n+'^H— > n + n+^K decay is 
not realized. As for the states ^'(21^-4/^,2/^)^1 which are dominated by the component '0"2~g^ (a, b) , they may 
be thought of as exhibiting some nn correlations in the internal region v — k < 5. 

In regard to the "democratic decay" of the ^H, it is highly improbable due to the fact that all limit eigen- 
values of the nucleus, namely Xf^, are different. Distinction of limit eigenvalues implies that there is no 
region where the Pauli-allowed states of the three-cluster system '^H coincide with hyperspherical harmon- 
ics, which reproduce the "democratic decay". Finally, Figure [3] demonstrates the coefficients C'(^i,_4^ 2/i)fc 
appearing in the expansion of the wave function T^(£;)(a, b) with energy E — 10.67 MeV above the 
^H^'^H-f n -I- n decay threshold in the basis of the Pauli-allowed states ^'(2i/-4/j. 2/1)1- • These coefficients have 
been obtained by diagonalization of the Hamiltonian, in which only the operator of the kinetic energy of the 
relative motion (with its exchange part) is retained. First five families were taken into complete account, 
families with fc > 5 were disregarded. Oscillator length ro was chosen so to reproduce root-mean-square 
radius of the "^H cluster. In fact, the eigenfunctions of the kinetic energy operator modified by the Pauli 
principle demonstrate the manner in which the antisymmetrization influences the expansion coefficients 
C(2u-Afj. 2fi)k' Figure [3] sustains our conclusion that families characterized by odd values of k contribute 
much more significantly to the wave function of the than families with even k do. An hierarchy among 
the coefficients of different branches of a given family is established by the behavior of their eigenvalues, 
namely, by the magnitude of attraction or repulsion in the corresponding SU(3) branches. By this reason, 
the coefficients of the branches (2i^, 0)/c=i take the lead, followed by the {2i',0)k=3 and {2iy — 4,2)^=3. Ne- 
glecting families of the states with k — 2 and k = A does not change considerably the wave function of the 
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Fig. 3. Coefficients CfL , „ , of the expansion of the wave function in the SU(3)-basis at i? = 10.67 MeV. Curves: 1 
- (2i.,0)fc=i; 2 - {2u,0)k^3; 3 - (2v-4,2)fe^3; J, - {2i/, 0)^=5; 5 - (2!/ - 4, 2)^=5; 6 - (2/^-8,4)^=5; 7 - {2u,0)k^2\ 8 - 
(2u-4,2)k=2; 9 - (2i.,0)fc=4; 10 - (2!.-4,2)fc=4; n - (2i/-8,4)fc=4- 

and the energy E. The latter takes such a large value mainly due to the fact that energy e of the "'H 
subsystem is high enough at given localization of the "'H. Of course, taking into account nucleon-nucleon 
interaction between cluster and a neutron would reduce this energy. Involvement of the "odd" families 
with higher k also leads to better description of the *H subsystem and, hence, to lowering the energy E. It 
should be also pointed out that the behavior of the expansion coefficients C^^u-a^ 2^)^ provides support for 
the predominance of the decay n+'^H^ n + n+'^H. 

6. Conclusion 

New approach to the problem of multichannel continuum spectrum of "^X = "^^^X + n + ^i systems [A < 6) 
is suggested based on the discrete representation of a complete basis of allowed states of the multiparticle 
harmonic oscillator, which was systematized with the help of the indices of SU(3) symmetry and defined 
in the Fock-Bargmann space. Proposed approach allows correct description of three-cluster systems both in 
the region of small intercluster distances, where the Pauli exclusion principle is of first importance, and in 
the asymptotic region where the scattering matrix elements are produced. 

Careful analysis of the structure of the eigenfunctions and behavior of the eigenvalues of the three-cluster 
norm kernel has been performed for the first time. A set of linear algebraic equations which generates the 
Pauli-allowed states for different three-cluster systems composed of an s-cluster and two neutrons is written 
in a general form. The latter depends on the single parameter which can be easily found for any three-cluster 
system. 

In the Fock-Bargmann space the Pauli-allowed states of a three-cluster system are shown to be superpo- 
sitions of the hypergeometric functions, with the expansion coefficients being orthogonal polynomials of a 
discrete variable. Eigenvalues of the three-cluster system are shown to tend to eigenvalues of a two-cluster 
subsystem with increasing the number of oscillator quanta i'. At the same time, corresponding eigenvectors 
take simple analytical form, while the expansion coefficients become the Kravchuk polynomials as the num- 
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ber of oscillator quanta increases. We suggest a way of resolving the problem of the SU(3) degeneracy of 
the Pauli-allowed states. A degree of the Kravchuk polynomial serves as an additional quantum number to 
label the states belonging to the same SU(3) irreducible representations. 

The Pauli-allowed states of a three-cluster system "^"^X+n+n can be arranged into branches and families, 
with all the states of a particular branch having common SU(3)-symmetry index /i, but differing in value of 
the first SU(3) index A. The eigenvalues belonging to a given branch tend to the same eigenvalue of the two- 
body subsystem Aj. ^ with the number of quanta increasing. The branches which share limit eigenvalues 
are combined in the family of the eigenstates, which thus is completely determined by the degree k of the 
corresponding Kravchuk polynomial. Each family of the Pauli-allowed states asymptotically corresponds to 
a certain binary decay channel of a three-cluster system into a two-cluster subsystem "^"^X occurring in a 
ground or an excited harmonic-oscillator state and a remaining neutron. Such asymptotic behavior gives an 
indication of possible decay channels of a three-cluster nucleus and allows us to specify the most important 
decay channels of the nucleus under consideration. 

The families which characterized by eigenvalues approaching limit values A^ from above are shown to 
be more favorable at small v, which can be considered as effective attraction of clusters due to the ex- 
change effects. It is precisely this interaction that causes the ^X system to decay via an intermediate stage, 
i.e.,^X ^^-1 X + n ^^"^ X + n + n. Alternative decay channel ^X ^^"^ X -l-^ n ^^"^ X -|- n -|- n is not 
realized. 

It is observed that distinction of limit eigenvalues of the three-cluster norm kernel results in the absence of 
such region where the Pauli-allowed states of the ^^^X+n+n system coincide with hyperspherical harmonics, 
which reproduce the "democratic decay". Hence, instead of hyperspherical functions the angular- momentum 
basis functions, which are labeled by the number of quanta, the angular momentum of the ^~^X subsystem 
and momentum of the relative motion of this subsystem and a remaining neutron, should be used in the 
limit V ^ k. As the result, asymptotically the three-cluster Schrodinger equation can be reduced to a two- 
body-like multichannel problem, while the asymptotic form of the expansion coefficients in the Z-basis is 
expressed in terms of the Hankel functions and the scattering S'-matrix elements. 

The validity of these conclusions was illustrated with the '^B.+n + n configuration of the nucleus. In 
particular, the structure of the Pauli-allowed states of the nucleus was revealed to correspond to the 
subsequent decay ^H— > n-|-^H— > n + n+^H. 



Appendix A. The Kravchuk polynomials 



The Kravchuk polynomial K!'^\m) of a discrete variable m, specified on the interval Q < m < v and 
orthogonal with weighting function and norm dk 



m!(i.-m)!' ~ \\ k\{y - k)\ 
can be defined via hypergeometric function: 



dk 



(-1)" R;;^ k)\{.-m)\ 
(v-m-k)\\J m\k\ f 



'2F1 (^-k, -m; u-m-k + 1; -^^ , 
if m < u — k 



= < 



^tZllv. \f^W^ ■ q'^P''-'^2F, (-, + k,-u + m;-u + m + k + 1; -|) , 

if m> v — k. 



Note that z/, m and k are the natural integers, < k < u. 
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Relation of orthogonality for the Kravchuk polynomials looks like: 



m=0 fe=0 

Symmetry relation can be written in the following form: 

Pv—m 



(_l)'»+fc/C^P)(TO) 



dk 



dv- 



The Kravchuk polynomial has the limit 



lim4^)(m) 



dk \r^k\^-n2'pqv 
where Hfe(a;) is the Hermitian polynomial. 



(m-'pv\ r {m-pyf^ 



(A.1) 
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